Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed February 14, 2013 (MN UT^ style file v2.2) 



o 



<: 
O 

^• 

Oh! 

O 



> 

(N 

O 

m 

O 

m 



Constraining the initial conditions of globular clusters using their 
radius distribution 



Poul E. R. Alexander^* and Mark Gieles^'^ 

^Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CBS OHA, UK 
^Department of Physics, University of Surrey, Guildford, GUI 7XH, UK 



Accepted 2013 February 09; Received 2013 February 08; in original form: 2013 January 21 



ABSTRACT 

Studies of extra-galactic globular clusters have shown that the peak size of the globular cluster 
(GC) radius distribution (RD) depends only weakly on galactic environment, and can be used 
as a standard ruler We model RDs of GC populations using a simple prescription for a Hubble 
time of relaxation driven evolution of cluster mass and radius, and explore the conditions 
under which the RD can be used as a standard ruler We consider a power-law cluster initial 
mass function (CIMF) with and without an exponential truncation, and focus in particular on 
a flat and a steep CIMF (power-law indices of and —2, respectively). For the initial half- 
mass radii at birth we adopt either Roche-lobe filling conditions ('filling', meaning that the 
ratio of half-mass to Jacobi radius is approximately rh/rj ~ 0.15) or strongly Roche-lobe 
under-filling conditions ('under-filling', implying that initially rh/rj ^ 0.15). Assuming a 
constant orbital velocity about the galaxy centre we find for a steep CIMF that the typical half- 
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light radius scales with galactocentric radius Rq as . This weak scaling is consistent with 
observations, but this scenario has the (well known) problem that too many low-mass clusters 
survive. A flat CIMF with 'filling' initial conditions results in the correct mass function at old 
ages, but with too many large (massive) clusters at large Rq. An 'under-filling' GC population 
with a flat CIMF also results in the correct mass function, and can also successfully reproduce 
the shape of the RD, with a peak size that is (almost) independent of Rq. In this case, the peak 
size depends (almost) only on the peak mass of the GC mass function. The (near) universality 
of the GC RD is therefore because of the (near) universality of the CIMF. There are some 
extended GCs in the outer halo of the Milky Way that cannot be explained by this model. 
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1 INTRODUCTION et al. (1991) and McLaughlin (2000) find R°c^ and R°g'^, respec- 

tively. Similar and shallower correlations are found for clusters in 
It has Ions been known that the shape of globular cluster (GC) lu- , , ■ nO 38 /■ i i i ■ i 

° f o external galaxies: rh oc Rq (in the Sombrero galaxy; Spitler 

minositv function IS relatively insensitive to galactic environment. . , ~„r,^^ r>o i7 ,r ^, ■ ■ r xt.^/-. cno 

et al. 2006), rh oc Rq (for the inner regions of NGC 5128; 
Accordingly, the peak luminosity of Mv — —7 is often used as „, „ ,, nO ii /■ i r ■ i 

, , „ , „ J „ , . Gomez & Woodley 2007), rh oc Rq (in a sample of six el- 

a standard candle (Shapley & Sawyer 1927; Reikuba 2012). It is „ , „„,..,.,. , . 

lipticals; Harris 2009). It clusters fill their Jacobi radius rj, their 
however not known whether this typical luminosity IS the result of ... , ■ jj- j ,i • i r r-. 

■' radii scale with their mass M and their orbital frequency \l as 

dynamical evolution (Fall & Zhang 2001) or the outcome of the » j-i/3r-.-2/3 tt c ■ » i r j i,-. i i 

■' _ ^ b > M ' S2 ' . Hence, for a constant M and orbital velocity Vg we 

cluster formation process (see Brodie & Strader 2006; Portegies 2/3 „ 

r, », „ D ^- 1 r^n^n c ■ , cxpcct OC RJ, , i.c. a stronger dependence on Rq than ob- 

Zwart, McMillan, & Gieles 2010, for reviews). f ni g ' b v 

served 

Jordan et al. (2005) showed that the half-light radii (rhi) of 
GCs are also sufficiently independent of galactic environment to be The objective of this work is to understand how the shape of 

used as a standard ruler Typical radii of ~ 3 pc are found for the cluster radius distribution (RD) depends on the cluster initial 
clusters in the Milky Way (Harris 1996, 2010 version) and other mass function (CIMF), the distribution of clusters within the galaxy 
galaxies (Kundu & Whitmore 2001; Jordan et al. 2005; Peng et al. and the initial cluster sizes. This will allow us to use the RD to put 
2011). Within galaxies the typical radii depend only weakly on additional constraints on the initial conditions of GCs. 
galactocentric radius Rq. For the Milky Way GCs van den Bergh 

The structure of this paper is as follows: In Section 2 we dis- 
cuss the RD arising for a synthetic sample of GCs at one Rq. In 
Section 3 we consider an entire population of synthetic clusters 
* e-mail: pera@ast.cam.ac.uk and compare our results with the Milky Way GC system. We finally 
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compare our results to previous studies and provide a summary in 
Section 4. 



2 THE RADIUS DISTRIBUTION OF CLUSTERS IN A 
SINGLE TIDAL ENVIRONMENT 

We start by considering tlie radius distribution (RD), i.e. tlie number 
clusters with a half-mass radius between rh and rh + drn, of clus- 
ters evolving in a fixed tidal environment (for example, evolving at 
a constant Rg within a galaxy halo). For this situation the RD can 
be related to the cluster initial mass function (CIMF, i.e. the num- 
ber of clusters with initial mass in the interval Mq and Mq + dMo) 
as 



dN 
dMo 



dM 



"3 — (-Kg) 
orb 



(1) 



Here \dMo/dM\ and \dM/dr\{RG)\ depend on the mass 
and radius evolution, respectively. These will be discussed in the 
following sections. 

2.1 Initial mass function and mass evolution 

In this section we assume that the CIMF is a power-law distri- 
bution: dN/dMo (X Mi" in the range Miow ^ Mo < Mup. 
In Section 3 we additionally consider Schechter (1976) distribu- 
tions. We aim to reproduce the GC population of a Milky Way like 
galaxy, and therefore adopt a constant Vg = 220 km s~^ such that 
Q oc Rq^- Throughout this work we assume that clusters evolve 
with a constant mass-loss rate that depends on the galactocentric 
radius as MRq = -20 MoMyr'^kpc (Gieles, Heggie, & Zhao 
201 1). The mass then evolves as 



M = Mo - A, with A = -Mt cc R(; 



(2) 



Here, t is the age of the cluster, for which we adopt 13 Gyr through- 
out this work. We then have \dMo/dM\ = 1. The total lifetime, 
or dissolution time, is tais = — Mo/M. For the final term of equa- 
tion (1), we need to know how the radii of clusters depend on M, 
for which we compare two different scenarios. In Section 2.2, we 
consider Roche-lobe filling clusters, whereby the ratio of half-mass 
to Jacobi radius {TZ = r-h/rj) is constant and approximately the 
Henon (1961) value, TZ = 0.145. In Section 2.3, we consider ini- 
tially Roche-lobe under-filling clusters, where initially TZ = Q and 
TZ increases during the evolution. 

2.2 Roclie-lobe filling clusters 

We first consider the scenario whereby all clusters fill their Roche- 
lobe from birth. Such clusters will evolve at a constant density 
whilst losing mass (Henon 1961) and so 



M = Art with A = oc Ra\ 



(3) 



where G is the gravitational constant. From this we find 
\dM/dr-h{RG)\ oc r^Ro^, which we combine with equation (1) 
to derive the RD 
dN , 



^{Ra)^{Ari + A) r^R^^ 



,.3«+2^-2(a+l) ^ A < Arg 



(4) 



For a = the RD is a single power-law with slope +2, while 
for a 7^ it is a double power-law distribution: at small radii the 
shape of the RD is independent of the shape of the CIMF and only 
depends on how clusters lose mass. This is because the Globular 
Cluster Mass Function (GCMF) evolves to dN/dM ~ constant at 
low masses, regardless of the shape of the CIMF (Henon 1961; Fall 
& Zhang 2001). At large radii the shape of the RD depends on 
the index of the CIMF, for which we find two regimes: firstly, for 
a < —2/3 the RD is an decreasing function of rh and the peak 
of the distribution occurs where the two power-laws join, i.e. at 
fh ^ {A/Ay^'^ oc Rq ^, where we introduce the symbol fh to 
denote the mode of the RD (most probable radius). For a > —2/3 
the RD rises until the largest radius present, which is the radius of 
the most massive cluster (Mup) and thus fh oc Rq^. 

2.3 Roche-lobe under-filling clusters 

We next consider clusters that are bom Roche-lobe under-filling. 
We use the model of Gieles et al. (201 1) for the relaxation driven 
evolution of rh as a function of Mo and Rg, in which clusters are 
assumed to form with a half-mass density much greater than the 
tidal density. In the first (roughly) half of their life clusters expand 
and become increasingly Roche-lobe filling, while in the second 
half they contract at a constant density, and near final dissolution 
the radius evolves as in equation (3) from Section 2.2. 

The functional form for the evolution of rh is (equation (26) 
in Gieles et al. 2011) 



rh ■■ 



My/3^^_ 



M 



1 5/2 



Mo 



2/3 



(5) 



Note that A = A{Rg), M = M(Mo, A) and A = A(i?G,t). 
In the late stages of evolution (Af <g; Mo), this function con- 
verges to the Roche-lobe filling relation (equation 3). Meanwhile, 
for M ~ Mo we find rh oc M^^^^t'^^^, which follows from the 
fact that in the expansion phase clusters evolve toward a constant 
ratio of Trh/i with Trh oc (Mfh)^/'^ (assuming a constant mean 
stellar mass and Coulomb logarithm, Spitzer & Hart 1971). If there 
are clusters in both the expansion and contraction phases, the size 
of the largest cluster (rJJ""') is not that of the most massive clus- 
ter, as was the case for Roche-lobe filling clusters (Section 2.2). To 
find r™^'' in a population of clusters with different Mo at a given 
Rg and t, we numerically solve drh/dMo = (equation 5), and 
find Mo(rh = r™ ^ 3.175A. We note that this is not the mass 
of the cluster that has just reached its maximum size, but instead 
a slightly more massive cluster that is still expanding (see Fig. 1). 
This is because the maximum rh is reached roughly half-way the 
total lifetime (Gieles et al. 201 1), so the cluster that has reached its 
maximum size has an initial mass of Mo ~ 2A. From Fig. 1 we 
also see that clusters can be placed in three categories, which we 
base on the their Mo : 

(i) Mo 5j 2A - clusters have past their individual maximum 
rh and are now approximately Roche-lobe filling (dashed/green 
tracks). 

(ii) 2A < Mo ^ 3.175A - clusters that are expanding toward 
their individual maximum rh, and show a positive correlation be- 
tween mass and r\ (solid/magenta tracks). 

(iii) Mo > 3. 175 A - clusters that are expanding, and have a 
negative correlation between M and rh (dot-dashed/cyan tracks). 

We cannot derive an easy analytical expression for the RD, 
because we can not invert the relation rh (Mo) (equation 5). We are 
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Figure 1. Radius evolution of clusters of different masses (10-* Mq < 
M < 10'' Mq) at i?G = 8.5kpc, Vq = 220kms-\ as a function of 
fractional dissolution time (left panel) and physical time (right panel). The 
thick/orange line represents the r\ of clusters after a Hubble time, and is 
curved in the left panel on account of the longer tjis of more massive clus- 
ters. The green/dashed, magenta/solid and cyan/dot-dash lines represent the 
evolutionary tracks of individual (sample) clusters in states (i), (ii) and (iii), 
respectively. These tracks are thicker for evolution that has occurred during 
a Hubble time, that is, on the left of the orange line. 



however able to give the behaviour of the RD at the extremes of 
rii. For expanding clusters, we use rh (x A/q ^^'^ and AI ~ Mg to 
find \dM/dri^\ oc r^'^, while contracting clusters approximately 
follow the same relationships as in Section 2.2. We substitute these 
expressions into equation (1) and find 

-(2+q) 

f4) / . , N , ^ TT,;,^ W 



- — (-Rg) oc 



-(3a 



,0 ==;rh <r£^ 
,rh(Mup) < rh < ri. 

Here we assumed that clusters with M = i\fiow have dissolved 
already and that rh shrinks to near dissolution. These two parts 
of the RD both occupy the region < rh < r™""" oc iij/^. The first 
approximation applies to the clusters whose rh has approximately 
reached the Roche-lobe filling condition {TZ ~ 0.145), i.e. A/o ^ 
3.175A, while the second approximation holds for the expanding 
clusters (Afo ^ 3.175A). For a — —2 we find that both extremes 
of the RD behave as rf-^. For very high Mup (where rh(Afup) ~ 0), 
this RD has the same shape as the RD for Roche-lobe filling clusters 
with a = (Section 2.2), but has a smaller r™'"'. 

As for the Roche-lobe filling case, the mode of the RD is 
either at the radius of the most massive cluster or at the radius 
of the cluster with A^o — few x A, depending on whether the 
RD of expanding clusters is rising or falling. The difference with 
the Roche-lobe filling case is that rh(Afup) is smaller than r™''". 
For a > —4/3 the RD of the Roche-lobe filling clusters is ris- 
ing and fh = r^'''' oc For a < -4/3 this part of the RD 
is a declining function, and if the most massive cluster is still ex- 
panding the RD peaks at the radius of the most massive cluster: 
rh = rh(Afup) < r™''" = constant (Fig. 2). This is the most 
straightforward way to explain a near universal radius scale for 
GCs. The condition that the most massive clusters have not yet ex- 
panded until their rh — TZrj is not satisfied in strong tidal fields, 
and will be discussed further in Section 3.4. 

The RD is computed for Roche-lobe under-filling clusters and 
for three values of a using a Monte Carlo approach (Fig. 2). The 
power-law behaviours at the expanding and contracting extremes 
are shown as dashed/black and dotted/red lines, respectively. The 
power-laws do not follow the Monte Carlo data where rh — r™^", 
as here Mq ~ 3. 175 A and mass loss is neither dominant nor neg- 
ligible. 

Now we have explored the different scenario for clusters at a 
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Figure 2. Monte Carlo simulations of clusters with different CIMF slopes 
(a) at Rq = 8.5 kpc with Vq = 220kms~^. The colour coding is the 
same as in Fig. 1 (green = stage (i) clusters, magenta = stage (ii) clusters, 
cyan = stage (iii) clusters). In the left panels, clusters in the three phases 
of evolution are shown separately, while they are added in the right panels. 
The dashed/black line represents the expected slope of expanding clusters 
while the dotted/red line corresponds to the expected slopes of approxi- 
mately Roche-lobe filling (contracting) clusters (equation 6). The sample 
size is 10^ clusters and clusters are evolved for 13 Gyr 



single Rq we will consider GC populations in a Milky Way type 
galaxy in the next section. 



3 RADIUS DISTRIBUTION OF A CLUSTER SYSTEM 

In this section we consider the RD of clusters in a Galactic halo, i.e. 
with a range of Rq . The RD can be found be multiplying the RD at 
a single Rq with the number of clusters at each Rq {dN/dRa ~ 
4TvRQn{RG), where n(i?G) is the number density of clusters in 



the halo) and integrating over all 7?g 



diV 
drh 



diV diV 
d^^'^^^d^R^ 



diiG. 



(7) 



The analytic evaluation of this integration gives results only for a 
limited number of cases, so we proceed by generating RDs in a 
Monte Carlo fashion. We therefore perform a number of simula- 
tions, for which we adopt a Hemquist (1990) model (i.e., a dou- 
ble power-law density profile with an central cusp and an r"* 
outer halo) with a scale radius a = 5 kpc and a truncation at Rq ~ 
100 kpc for the distribution of clusters in the Galaxy {n{RG)), 
which is a reasonable approximation for GCs in the Milky Way 
(McLaughlin 2000). For the CIMF we use a Schechter (1976) func- 
tion: dN/dMo oc Afo" exp(-Afo/A/,), with M, = 10^ Mq and 
Afiow ~ 100 Mq. For the slope we consider a = —2, i.e. what 
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is found for young clusters (Zhang & Fall 1999; Portegies Zwart 
et al. 2010) and a = 0, which describes the MF of old clusters in 
the Milky Way and early type galaxies in the Virgo Cluster (Jordan 
et al. 2007). A CIMF with a = is invariant under mass-loss 
(Jordan et al. 2007) so at an age of 13 Gyr the GCMF has the same 
shape as the CIMF. The resulting RDs and the fh {Rg) relations for 
these four cases are shown in Fig. 3 and are discussed below. We 
compare all these simulations to the rh values of Milky Way GCs 
compiled from the (Harris 1996, 2010 version) catalogue. Values 
for rh are derived by multiplying the observed half-light radii by 
4/3 to correct for the effect of projection. 

3.1 Scenario A: Roche-lobe filling clusters and cx.— —2 

We first consider Roche-lobe filling clusters and a CIMF with slope 
a = — 2. We recall from Section 2.2 that for Roche-lobe filling 
clusters there is a critical value for the CIMF slope of a = —2/3: a 
steeper CIMF gives rise to a double power-law RD, rising at small 
radii and falling at large radii. The total RD is the sum of many such 
functions, and therefore has a (logarithmic) slope of +2 at small 
radii. At large radii the integration of equation (7) gives a slope of 
3a -t- 3.5 = —2.5. The typical radius scales as fh oc R]!'\ We 
recall from Section 2.2 that this is due to the scaling of the break in 
the double power-law RD with _Rg. This scaling is close to what is 
found for Milky Way GCs, but in this scenario there are too many 
large clusters (see panel A, Fig. 3). In addition, too many low-mass 
clusters survive at large Rq in this scenario (Vesperini 2001 ; Gieles 
& Baumgardt 2008). 

3.2 Scenario B: Roche-lobe filling clusters and a = 

For Roche-lobe filling clusters and a = the RD at a single _Rg 
is simply a +2 power-law. In this case the integration of equa- 
tion (7) can be performed, with the result that at large radii we 
find dA'"/drh on r^^ . Because of this slow decline, this scenario 
produces even more large clusters than scenario A. However, ow- 
ing to the invariance of this CIMF under the influence of mass-loss, 
this scenario does predict the correct GCMF (meaning that the peak 
mass, M, is independent of i?G, Vesperini 2000). Because fh scales 
with the Jacobi radius of the most massive cluster (Section 2.2) at 
any given Ra, we find fh oc R^^ (see panel B, Fig. 3), which is a 
stronger correlation than found for Galactic clusters. We note that 
this is the only scenario our of the four we considered for which the 
scaUng is foimd. 

3.3 Scenario C: Roche-lobe under-filling and <x— —2 

For under-filling clusters with a steep CIMF (a = —2) both the 
expanding and contracting clusters have a RD that behaves as a 
-1-2 power-law (just as in scenario B). However, the total RD is 

different because the typical radius scales with fh oc r™^^ and 
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therefore fh oc Rq (Section 2.3). Accordingly, there are fewer 
large clusters in this scenario, because high mass clusters are still 
expanding to their tidal boundary and are smaller than clusters of 
equivalent mass in scenario B. 

3.4 Scenario D: Roche-lobe under-filling and a — 

In the final scenario (under-filling clusters with a > —4/3) the 
sample population contains a fraction of expanding clusters, which 



is larger for higher values of a (Fig. 1). The radii of expanding clus- 
ters are independent of their local tidal field and also independent 
of their initial radii (Henon 1965; Gieles et al. 2010). This is the 
case for _Rg ^ 5 kpc (panel D, Fig. 3). The typical n is then sim- 
ply the radius of the cluster with the typical mass. For _Rg < 5 kpc, 
the majority of clusters have expanded to become Roche-lobe fill- 
ing and are therefore described by scenario B with a typical radius 
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scaling as f h oc Rq . The number of clusters in this regime de- 
pends on the slope of the CIMF and the Hemquist scale radius. 
This scenario is the only scenario that produces a clear peak in the 
radius distribution below 5 pc, and the only scenario to produce a 
rapid decline in the RD at high rh consistent with Milky Way GCs. 



3.5 Comparison to Milky Way GCs 

In the right panels of Fig. 3 we show the RD and the individual rh 

against Rg values of the MW GCs. The width of the Milky Way 
GC RD is narrow, with radii rarely in excess of 20 pc (90% of the 
GCs have radii < 10 pc, Harris 1996, 2010 version). Comparing 
this to the models we see that in Scenarios A, B there are too many 
clusters with rh > 20 pc (see Fig. 3), especially in scenario B, 
where we find clusters with size up to rh ^ 100 pc. Scenarios C 
and D both show a narrower range of rh, in better agreement with 
that found for Milky Way GCs, with scenario D providing a very 
comparable shape overall. 

Several MW GCs at high Rq are larger than those predicted. 
These outlying clusters cannot be explained by our simple model 
for cluster evolution, and so may result from effects neglected by 
this study (e.g., these clusters may be accreted, or perturbed by the 
tidal fields of external galaxies). The scaling of f h with _Rg of the 
Milky Way GCs shows considerable scatter, but is approximately 
rh oc -Rg*, similar to the scalings of scenarios A and C and shal- 
lower than scenario B. However, the majority of Milky Way clus- 
ters are found within the innermost ~ 10 kpc. In this regime the 
radii in scenario D are similar to those in C. It is therefore difficult 
to satisfactorily distinguish which scenario optimally represents the 
rh scaling of Milky Way clusters, but combined with the superior 
match of the RD shape we conclude that scenario D is preferred. 
The rh(J?G) relationship of scenario D is also similar to that de- 
scribed for M87 in Webb et al. (2012). We therefore conclude that 
the good match of the shape of the RD in scenario D suggests that 
most clusters form Roche-lobe under-filling, and with a relatively 
flat CIMF, similar to the present day GCMF. 



4 SUMMARY AND DISCUSSION 

We have investigated the radius distribution (RD) of globular clus- 
ter (GC) populations in Milky Way-type haloes for a variety of ini- 
tial conditions for the masses and radii. We find that Roche-lobe 
filUng initial conditions results in too many large clusters (> 20 pc) 
at large Rg, for both steep and shallow cluster initial mass fimction 
(CIMFs). Models in which clusters are initially Roche-lobe under- 

1/3 

filling and form with a steep CIMF show a scaling of f h oc Rq , 
roughly consistent with Milky Way GCs, but the shape of the syn- 
thetic RD is wider than the observed one. The mode of the half- 
mass RD of a GC population is independent of Galactic environ- 
ment only if clusters form Roche-lobe under-filUng, and with a rel- 
atively flat CIMF. 

Alternative models to explain the characteristic GC radius of 
a few pc exist (e.g. Jordan et al. 2005; Harris et al. 2010). In these 
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Figure 3. Monte Carlo samples of N = 10^ (surviving) star clusters at t = 13 Gyr within an isothermal galactic halo. Clusters are simulated with a Hemquist 
(1990) distiibution of galactocentric radii and a Schechter (1976) CIMF (see text for parameters). In the top row cyan, magenta and green regions correspond 
to clusters undergoing the three evolutionary stages outhned in Section 2, whilst in the bottom row individual sample clusters are denoted by green diamonds 
(stage (i) clusters), magenta squares (stage (ii) clusters), or cyan dots (stage (iii) clusters). The left hand panels represent scenarios A, B, C, and D as labeled, 
while the rightmost figure corresponds to the present day Milky Way GC population compiled from Harris (1996, 2010 version). The red solid line in the 

1/3 2/3 

bottom row is the median rjj as a function of Rq, while the black (dotted) lines represent rjj = const, rjj oc Rq , and rjj oc R/ as appropriate. 



models the present day size is the resuU of the star formation pro- 
cess or the gas expulsion scenario, respectively. The subsequent 
Hubble time of evolution is assumed not to significantly change the 
radii of GCs. In our model scenario D, however, the typical GC size 
is the natural outcome of a Hubble time of dynamical relaxation 
and is not sensitive to the exact initial conditions or early evolution 
since the relaxation driven expansion erases the initial properties of 
the clusters (in terms of rh, Henon 1965). 

In addition to these works, Madrid et al. (2012) show (via A'^- 
body studies) that the radius is set "by clusters that are originally 
massive enough to survive mass loss due to tidal stripping". Thus, 
the present day radii in their model depends only on clusters' abil- 
ity to survive in the Galactic tidal field. They consider clusters at 
different Rq, all with the same mass and rh. Their final rh vs. Rq 
relation (their Fig. 3) is very similar to our rh(-RG) relation found 
in our scenario D. The radii of their clusters at large Rq are larger, 
because they consider clusters with lower masses which have ex- 
panded more (Henon 1965). 

In contrast to our results, Shin et al. (2013) find from a large 
suite of Fokker-Planck simulations that models in which GCs are 
on average larger than at the present fit the data best. A fraction 
of their clusters are Roche-lobe overflowing at birth. Because their 
result is based on fits against the entire GC population, exclud- 
ing accreted clusters, their fit results are mainly driven by clusters 
in the inner parts of the Milky Way. Also in our scenario D most 
clusters at Rq < 5 kpc have filled their Roche-lobes and it is hard 
to distinguish scenarios there. Analyses of the orbits and radii of 
GCs suggests that a large fraction of them are Roche-lobe filling at 
the present day (e.g., Baumgardt et al. 2010; Ernst & Just 2013), 
supporting our conclusion that most clusters started more compact. 
Shin et al. (2013) find that the CIMF must have been flat, in agree- 
ment with what we find. The flat CIMF needed to get the shape of 
the RD right is in stark contrast with the observed power-law mass 
function of young clusters (Portegies Zwart et al. 2010) and is yet 
another argument that dynamical evolution cannot be responsible 



for the difference (see also Vesperini 2001; Gieles & Baumgardt 
2008). 

For the models described in this study we assumed the in situ 
formation of the GC population. However, some of the outer halo 
clusters, many of which are the large outliers in the RD were most 
likely from dwarf satellites. A connection between dwarf galaxies 
and extended clusters has previously been suggested (Elmegreen 
2008; Da Costa et al. 2009; Mackey et al. 2010). We have further- 
more made the assumption of a constant mass-loss rate, and we 
omitted several physical effects such as the delayed escape of stars 
from the Roche-lobe (Fukushige & Heggie 2000) and stellar evo- 
lution. A fast cluster evolution code is being developed (Alexander 
& Gieles 2012) with the goal to include all these effects to be able 
to model the RD, GCMF and the Rq distribution in the future. 
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